Constructing local composite operators for glueball states 
from a confining Gribov propagator 



M. A. L. Capri a *, A. J. G6mez bt , M. S. Guimaraes 6 *, V. E. R. Lemes 6§ 

S. P. Sorella 6 1 , D. G. Tedesco b ** 

a UFRRJ — Universidade Federal Rural do Rio de Janeiro 
<^> Departamento de Fisica — Grupo de Fisica Teorica e Matemdtica Fisica 

BR 465-07, 23890-971, Seropedica, RJ, Brasil. 
^ b UERJ — Universidade do Estado do Rio de Janeiro 

Q Institute de Fisica — Departamento de Fisica Teorica 

Rua Sao Francisco Xavier 524, 20550-013 Maracana, Rio de Janeiro, RJ, Brasil. 



Oh 



Abstract 



The construction of BRST invariant local operators with the quantum numbers of the lightest glueball 
states, J PC = ++ , 2 ++ ,0 h , is worked out by making use of an Euclidean confining renormalizable 
gauge theory. The correlation functions of these operators are evaluated by employing a confining 
gluon propagator of the Gribov type and shown to display a spectral representation with positive 
spectral densities. An attempt to provide a first qualitative analysis of the ratios of the masses of the 
^2 lightest glueballs is also discussed 
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1 Introduction 



In a previous work [1], we have outlined the construction of an Euclidean nonabelian gauge model, called 
the replica model, in order to study aspects of the gluon confinement. As the name let it understand, 
the model relies on the introduction of a replica of the usual Faddeev-Popov action in the Landau gauge, 
where the gauge fields are coupled through a soft term, i.e. through a term which is quadratic in the 
fields. The model turns out to have many aspects in common with the Gribov-Zwanziger theory [2j[3lH], 
a feature which enables us to employ it in order to improve our current understanding of the role played 
by the Gribov horizon as well as of its consequences on the analytic properties of the correlation functions 
of gauge invariant local composite operators. 

As reported in [1], the replica model enjoys the following properties 

• the gluon propagator turns out to be a confining propagator of the Gribov type, namely 

W*)4(-*)> = (<W - kj ^) , (i) 

where the mass parameter $ 2 plays a role akin to that of the Gribov mass parameter 7 2 of the 
Gribov-Zwanziger theor}|^] 

• as the Gribov-Zwanziger theory, the replica model also displays a soft breaking of the BRST sym- 
metry 

• it exhibits a nice interpretation in terms of i-particles [6], i.e. in terms of the two unphysical modes 
with complex conjugate imaginary masses ±i\/2i? corresponding to a confining Gribov propagator 

k 2 I ( 1 1 

fc 4 + 2tf 4 ~ 2 U 2 -n/2tf 2 + k 2 + iy/2d 2 

As shown in |6j, the introduction of the i-particles provides a powerful tool in order to construct 
examples of local composite operators whose correlation functions do have good analyticity proper- 
ties, as expressed by the Kallen-Lehmann spectral representation, and this in spite of the fact that 
they are evaluated by using a Gribov type propagator, eq. ([!]), 

• analogously to the case of the Gribov-Zwanziger theory [8], the replica model admits a refined 
version [I], giving rise to a decoupling type gluon propagator which does not vanish at the origin 
in momentum space, k 2 = 0, i.e. 

W*>4(-*» - y 2 ^ (V - ffi ■ (3) 

In the present paper we pursue the analysis of the replica model. We prove that the model is renor- 
malizable to all orders and we discuss the construction of a set of local BRST invariant composite 
operators whose correlation functions exhibit a spectral representation. We focus, in particular, on the 
BRST invariant operators having the quantum numbers of the three lightest glueball states^J namely, 
J PC = ++ , 2 ++ , ^ . The corresponding spectral densities are worked out at one- loop order and shown 

See ref. [5] for a general introduction to the Gribov-Zwanziger theory and its refined version. 
2 See ref. [5] for a recent review on glueballs. 
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to be positive. A first qualitative analysis of the ratios of the masses m^ ++ , m^ ++ , ?tIq-+ is also presented 
by relying on a phenomenological SVZ-type sum rules0 [TO H21 E9 E] from which an expression for the 
glueball masses is established through the use of the Borel transformation. 

The work is organized as follows. In Sect. 2 the main features of the replica model are briefly reviewed. In 
Sect. 3 we derive the set of Ward identities fulfilled by the model and we establish its renormalizability to 
all orders by using the algebraic renormalization |15j . These results are proven to extend to the refined 
version of the model. Also, a simple argument accounting for the non-renormalization of the mass param- 
eter # 2 appearing in the Gribov propagator, eq.Q, is provided. Sect. 4. deals with the construction of 
local BRST invariant operators with the quantum numbers of the lightest glueball states ++ , 2 ++ , . 
The relationship of these operators with the z-particles is discussed and the Kallen-Lehmann spectral 
representations for the corresponding correlation functions is established at one-loop order. A qualitative 
analysis of the ratios of the masses m^, m Q-+ is reported in Sect. 5. Finally, Sect. 6. gathers our 

conclusion. 



2 A short survey on the replica model 

One starts by considering the Faddeev-Popov action in the Landau gauge 

S FP = J d 4 x (\F^(A)F^iA) + ib a dpA% + ^Df(A)A , (4) 

where b a is the Lagrange multiplier enforcing the Landau gauge condition, d^A^ = 0, (c a ,c a ) stand for 
the Faddeev-Popov ghosts, Ffi v (A) is the field strength 

F; V {A) = d,Al - d„A° + 9 r**A\A% , (5) 

and D^ b (A) denotes the covariant derivative 

DfWc^d^-gr^A^. (6) 

The action of the replica model is obtained as follows: 

• one first considers a replica of the Faddeev-Popov action, eq.Q, by introducing a set of mirror 
fields (U®,b a ,u) a ,uj a ) as well as a mirror Faddeev-Popov action: 

Smfp = J d 4 x (\u^(U)U^(U) + ib a d,U; + O^DfiU)^ , (7) 

where U^ U (U) is the field strength corresponding to the field 

U; U {U) = d^l - dJJl + gf abc U*U c „ , (8) 

while b a and (O a ,uj a ) stand for the mirror Lagrange multiplier and mirror Faddeev-Popov ghosts, 
respectively, and 

Df(U)u b = dpu* - gf abc U c ^ b . (9) 

3 See ref.[TD] for a pedagogical introduction to the SVZ sum rules and their applications to QCD. 
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the gauge field A a is softly coupled to its mirror field Z7? through the following mixed term 



S$ = iV2$ 2 J d A x A*U% , 



(10) 



where, as already mentioned, i9 2 is a mass parameter which plays a role analogous to that of the 
Gribov parameter 7 2 of the Gribov-Zwanziger action [21 [3111] • 

The replica model is thus specified by the following action 

Sreplica = $FP + SmFP + S$ 
r j4 



+ib a d„Al + c a d^Df{A)c b + fb a d„U* + uo a d^D ab {U)u l 



(11) 



Let us have a look at the propagators of the fields (A a , J7") 



i.e. 



«6 



(A«(k)Al(-k)) 
(U a Jk)U b u (-k)) = 6 ab 



(A a Jk)U b u (-k)) = 8 ab 



k 2 




k^kp 


fc 4 + 2tf 4 


k 2 


k 2 


V - 

V - 


k^k v 


fc 4 + 2tf 4 

-iV2# 2 
k A + 2tf 4 


k 2 

k^ky 
k 2 



(12) 

(13) 
(14) 



As one sees from expressions ( |12[ ), (13), (14), all propagators are of the confining Gribov type. As such, 
they correspond to unphysical excitations. Said otherwise, we cannot attach a particle interpretation to 
the propagators of the elementary fields (A^, UT). 

Another feature displayed by the model is that it contains a unique coupling constant g. Both fields 
A^ and U£ interact with the same coupling. We remark that the feature of having a unique coupling 
constant is protected by a powerful discrete mirror symmetry pj. It turns out in fact that the action 
Sreplica, eq.(ll), is left invariant by the following discrete transformations 



ACL _ 

u a - 
b a 
If 



id 



>b a 

> b a 
-oj a 

> c a 
-uj a 



(15) 



The mirror symmetry (15) means that the fields (A®, b a , c a , c a ) can be replaced by ([/", b a , uj a , oj a ), and 



vice- versa. As in the case of the Gribov-Zwanziger theory [5], the action (11) displays a softly broken 
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BRST invariance. It turns out in fact that the nilpotent BRST transformations 





— —n ab (A)r b 




T\ab fjj\ , ,b 


sc a 


= 9 -f abc c b c c , 
2 


sw a 


= 9 -f abc uj b u c , 
2 


sc a 


= ib a , 


sb a 


= o, 


soo a 


= ib a , 


sb a 


= o, 



leave the action S rep u ca invariant up to soft terms proportional to the parameter t? 2 , i.e. 

sS re pH ca — l? A^reafc > 

where A 6reafc is given by 

Aftrea* = -*v^ | d 4 x (u*D ab {A)c h + A»Df(U)u b ) . (18) 

Being of dimension two in the fields, Ab rea k is a soft breaking. Finally, as discussed in [I], the replica 
model reduces to the Faddeev-Popov action when -i? 2 = 0. 



2.1 Relationship with the i-particles 



In view of the construction of the BRST invariant local composite operators with the quantum numbers 
++ , 2 ++ , it is worth to spend a few words on the relationship between the replica model and the 
i-particles [6j, i.e. the pair of unphysical modes with complex masses ±iy / 2'i? 2 associated to a confining 
Gribov type propagator. To this end, let us consider the quadratic part of the action S rep u ca containing 
the two gauge fields (A®, {/"), namely 

s guad = J d*x f^(-d 2 )A- + l -u;{-d 2 )u; + iV2# 2 A«u«\ , (19) 

where we have already taken into account the Landau gauge conditions, d^A® = and d^U^ = 0. 
Expression (21) can be cast in diagonal form by introducing the two field combinations 



< = 71 1£% ~ K) ■ (2°) 



Therefore 



S quad = fd A x QA^-9 2 + i v / 2 1 9 2 )A« + ^(-a 2 -iv / 2^ 2 )r ? ^ 



(21) 
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which describes in fact the propagation of two unphysical modes with complex masses ±ii/2'& 2 



-k)) 
-k)) 
-k)) 



^((A-(k) + U-(k))(A b v (-k) + U b u ( 



((A«(k)-UZ(k))(Al(-k)-Ul 
((A a Jk) + U a Jk))(A b u (-k)-U b u 



-k))) 
-k))) 
-k))) 



tab 



1 



k 2 + iV2$ 2 



sab 



k 2 - iV2$ 2 



'flU 



'\1V 



k^k v 
kpk v 



. 



(22) 



These are precisely the i-particles corresponding to the Gribov propagators in eqs.(12), (13), (14). We 
see thus that the action (11) has a direct interpretation in terms of «-particles. As observed in [6], the 
advantage of introducing the fields (A° , rf^) relies on the fact that they turn out to be helpful in order to 
construct examples of local operators whose correlation functions exhibit the Kallen-Lehmann spectral 
representation. This feature stems from the observation that the momentum integrals corresponding to 
Feynman diagrams containing an equal number of propagators of the A-type and of the ?7-type can be 
cast in the form of a spectral representation, see [6] for a detailed discussion. Let us consider in fact the 
one-loop integral 

dS 1 (23) 



X(k 2 



(2tt) 4 ((k - p) 2 + iV2$ 2 ) (p 2 - iV2# 2 ) 



For the corresponding spectral representation one has [6] 



l(k 2 ) -1(0) 



drp(r)(- 

J2V21) 2 V 



+ k 2 



where the spectral density 



P(r) 



Vt 2 



8tf 4 



(24) 



(25) 



16-7T 2 T 

is positive in the range of integration^} This relevant property enables us to construct local operators 
with good analyticity properties. As an example, we quote the operator 



(26) 



also extensively investigated in [6], where its two-point correlation function has been shown to be cast in 
the form of a spectral representation with positive spectral function^] 



{0 Xv (k)Ox v (-k)) 
P(r) 



r + k 2 



12(N 2 - 1) 



Vt 2 - 8tf 4 (8# + r 2 ) 
32vr 2 r 



(27) 



3 Renormalizability of the replica model and non-renormalization prop- 
erties of the mass parameter $ 2 

The existence of a soft breaking of the BRST symmetry does not prevent us to establish a set of Slavnov- 
Taylor identities which are suitable for an all-order algebraic analysis of the renormalizability properties 



4 The subtraction of the factor 1(0) in eq.(24l is needed to account for the divergent character of expression ( |23[ ) in four 
dimensions. 

5 Also here, a suitable subtraction to get rid of ultraviolet divergences is needed, see [B] for details. 
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of the model. The usual way of proceeding is that of introducing an extended action which incorporates 
all local composite operators entering the soft breaking, by coupling them to a suitable set of external 
sources. The original action is thus recovered when the sources acquire a particular value, which we shall 
refer to as the physical value. The renormalizability of the extended action entails thus the renormaliz- 



ability of the starting action S re piica, eq.(ll). 



The use of the so-called algebraic renormalization is particularly suited for the study of a field the- 
ory whose starting classical action is characterized by a set of symmetries, as in the case of the replica 
model considered here. For such theories the algebraic renormalization provides a purely algebraic set 
up in order to establish an all orders proof of the renormalizability in a regularization independent way. 
Essentially, the algebraic renormalization combines in a powerful way the locality properties of the per- 
turbative series, as expressed by the power-counting, with the set of Ward identities following from the 
symmetry content of the starting action. More precisely, from the power-counting, one knows that the 
ultraviolet divergences which originate through radiative corrections are local polynomials in the fields 
and their derivatives whose dimensions are bounded by four. These local polynomials are thus completely 
characterized by the requirement that they have to be compatible with the set of Ward identities corre- 
sponding to the symmetries of the starting action. In much the same way, the power-counting ensures 
that potential anomalies are also associated to local polynomials in the fields and their derivatives which 
are constrained by the Wess-Zumino consistency conditions stemming from the algebraic relations among 
the various functional operators associated to the Ward identities. In summary, both invariant countert- 
erms and potential anomalies can be fully characterized in a purely algebraic way as the most general 
solution in the space of the local field polynomials of the Ward identities which express the symmetry 
content of the model. This program goes under the name of algebraic renormalization |15j . 



3.1 Identifying the classical extended action 



In order to write down the Slavnov- Taylor Ward identities, we follow the general set up of the algebraic 
renormalization |15j and start by introducing a term depending on BRST invariant external sources 
coupled to the nonlinear BRST transformations of the fields appearing in eq.(16), namely 



^external 



d 4 x 
d 4 x 



n a JsA°) + L a {sc a ) + WJsUZ) + L a {suj a ) 



-Q a u D ab (A)c b + ^-f abc L a c b c c - O" D ab (U)u b + - f abc L a u b u' 



(28) 



Notice that the mirror symmetry (15) can be immediately extended to the BRST sources by requiring 
that 



5 V 



L a -> L a , 

We proceed now by introducing a BRST doublet 



sK 



ab 



J 



ab 



O" — \ Q a 



L a -> L a . 

of external souce^] 



sJ ab = . 



(29) 



(30) 



6 See, for example, the case of the Slavnov- Taylor identities derived in the Gribov-Zwanziger theory [5l IHl 116] . 
7 As we shall see later, the introduction of these sources with two color indices (a, 6) is useful in order to forbid the 
appearance of counterterms of the kind (A^A^ + U^U^). 
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and obtain the following extended BRST invariant term 

Sj = s j <fxK ab Apl = j d 4 x [j ab A°U b + K ab (Df(A)c c )U b + K ab A*D b °(U)oJ c ^ . (31) 



As one can easily check, the soft term S#, eq.(10), can be recovered from the invariant term Sj when the 
external sources J ab and K ab attain the following physical values 

J ab \ =iV2tf 2 6 ab , K ab \ , =0, Sj\.=St. (32) 

Ipnys 'phys 'phys v 7 



Also, the invariance under the mirror symmetry (15) is guaranteed by demanding that 



jab _^ jba ^ R ab _^ R ba _ (33) 

Finally, for renormalization purposes, an extra BRST exact term has to be introduced 

S* cxtra = \ S j (K ab J ab - gf*bc K ad K bd c c _ g fabc K da K db u c 

= C J d A x Q J ab J ab - gf abc J ad K bd c c - gf abc J da K db uj c - ^ f abc f cmn K ad ''K ba l c m c n 

_^fabcfcmn K da K db UJ rn UJ n\ _ ^ 



Here, ( is a constant parameter and, when the physical values of the sources (32 ) are taken, only the first 
term of the r.h.s. survives, giving rise to a constant vacuum term 

S extra \ phy = -((N 2 -l)^V, (35) 

where V is the 4-dimensional Euclidean volume. 

Thus, the starting extended action we shall consider reads 

E = SfP + -SmFP + S'external + Sj + S extva (36) 

3.2 Ward identities 



In oder to discuss the renormalizability of the extended action S, eq.(36), let us give the full set of Ward 
identities fulfilled by S. These are 

• the Landau gauge fixing equations [15] 

w =id ^ w= id ^- (37) 



the anti-ghosts Ward identities [15 



8 



the Slavnov- Taylor identities |15j 

<S(E) = , (39) 

with 



<?(y^ — / ( y 4 T ( I __ j I i ju a i_ ju a i_ j at > I 

v ' J \5A% <ffi£ 5U* <5J2» 5c a 5L a 5u a 5L a 5c a 5u a 5K ab J 



(40) 



the two rigid symmetries 



W(E) = gf*> l* x lA>£- + + c b f§ + c b ^ r + ill 



P6A?, 5b c 5c c 5c c * SQ^, 



6E r ,5S udH uSE ^ SE 



+Lb jL< +jM SJ* +KU SK*)=°' (41) 



These symmetries imply that the first color index of the external sources J ab and K ab can be 
contracted only with the fields of the ^4-family, i.e. (A't, b a , c a , c a , O", L a ), while the second index 
can be contracted only with the fields of the [/-family, (U®,b a ,u) a ,u} a ,£l^,L a ). They forbid, for 
example, the presence of counterterms like J^A^A^ and J ab U^U b . 

the two ghost Ward identities [TTJ, [15] 

£<*(£) = A» lass , g a (Z) = A a cl3SS , (43) 



where 



and 



o a = j d ^{^-^r b ^^ c +9f abc K bd ^), (44) 

A£ lass = J d 4 xf abc ^A c ^-L b c c -g 2 Cf amc f dbq K md K cb u q ^ , (45) 



^■[^-igf^&JL+gf^K*^-]. (iG) 



A^ lass = J d 4 X f ab ° (n^U* - L b 0J C - g ^f*mc f db qK dm K bc c , 



(47) 



Notice that the expressions A", , A"i in eq.(43) are linear in the quantum fields. As such, they 



represent classical breakings, not affected by quantum corrections |15| . 
the SL(2, M) symmetries 



2>(E) = l*x(f^-M£-)=0, (48) 
' V 5b a 5L a J v ' 

P(-) - /^f^-^I^Vo. (49) 
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3.3 Characterization of the most general invariant counterterm 



5b a c 


= o, 


5b a 


J* 


= o, 


(- 




= o, 


W a s c 




= o, 






= o, 






= o, 





In order to characterize the most general invariant counterterm which can be freely added to all orders in 
perturbation theory [J5], we perturb the classical action E by adding an integrated local polynomial S c 
of dimension bounded by four, and with vanishing ghost number. We demand thus that the perturbed 
action, (S + r/S c ), where rj is an expansion parameter, fulfills, to the first order in rj, the same Ward 
identities fulfilled by the classical action S. This requirement gives rise to the following constraints for 
the counterterm E r 

0, 

= 0, 
0, 

= o, 

(50) 

where 5s, 2?s and 2?s denote the linearized operators 

f 4 f SE S | <5 | <SS 5 <JS 5 5E 5 t SE 5 

s ~ J X \SA^59^ + SQ^SA^ + SU^W^ + Sn^6U^ + 5c"5L« + «^<5c" 

j i5S i a 5 - a 5 b 5 \ 

p s = d i x[c a - z — — -! — — , 52) 

7 V 5c a 5b a 5L a 5L a 5b a J v ' 

- f . ( _ SE .5E 5 . d£ 6 \ 

£> s = I d 4 x[uj a - z - - . (53) 

J \ Sufi <^ a <^^ a <5L a Sb a J v 7 

Following the general set up of the algebraic renormalization [15], it turns out that the most general 
invariant counterterm S c obeying the constraints ( 50 ) and the mirror symmetry is given by 

S c = «o J d 4 x (~F; v F; v + \u^U^\ + S s | d 4 x [ ai ((Q« + d^A^ + + d^ a )t£) 

+a 2 K ab A a ^U* + ^ (K afe J ab - gf<*bc K ad K bd c c _ g jabc R da R db y c) j _ (54) 

where ao, ai, 02,03 are free coefficients. 



3.4 Renormalization factors 



It remains now to show that the invariant counterterm ( 54 ) can be reabsorbed through a redefinition of 



the parameters, fields and sources of the classical starting action S, according to 



Z 



1/2 



(55) 



10 



where 



{^4, U, 6, b, c, u, c, oj} , 
{g,(,n,ti,L,L,J,K}, 



so that 



£(00, $o) = *) + »/S c (^, + 0(r? 2 ) 
By direct inspection, the renormalization factors are found to be 



Z 



1/2 



z, 



-1/2 



z 



-1/2 



z 



1/2 
A ' 



„l/2 _ yl/2 _ ryl/2 _ „l/2 _ r? _ ry 

— &c. — — ~ — "C 



Z -l/2 = zr l/2 
L L 



7 -l/2 7 -l/4 
An A » 



Zr = ZjZg 1 l 2 z)( 4: , 



with 



Z 



1/2 

Zj 



1 + e 



a + 2ai 



1 



a 



2 ' 

1 - e (a - a 2 ) , 

1 + e (2a - 2a 2 + a 3 ) . 



(56) 
(57) 



(58) 



(59) 



Equations (57), (58), (59) show that the counterterm S c can be reabsorbed by means of a redefinition of 



the fields, sources and parameters of the starting action E, establishing thus the renormalizability of the 
replica model. 



3.5 Non-renormalization properties of the mass parameter $ 2 



In this section we shall show, by graphical arguments, that the term 



M ' I phys 



does not renormalize, namely 



which implies that the mass parameter $ 2 enjoys the following non-renormalization properties 



Z,Q2 — Z A 1 , 



(60) 
(61) 

(62) 



and <i2 = —1a\ in eqs.(59), meaning that its renormalization factor Z$i is not an independent quantity. 



Let us first pay attention to the two-point 1PI Green's function of the gauge field A a ^ at one-loop ordeiQ 
In this case, a counterterm like ^A^A^ could arise, for example, from a tadpole diagram. This diagram 



is related to the following momentum integral 



d d q 



(2vr) d g 4 + 2tf 4 



(63) 



Due to the mirror symmetry, eq.(15l, equivalent results can be obtained for the replica field [/". 
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where d = 4 — e and where use of the dimensional regularization is understood. Noticing that the Gribov 
propagator displays the following property 

q 2 1_ 2^ 

" ^2 



q A + 2tf 4 

it follows that the integral I can be rewritten as 



q 2 q 2 (q 4 + 2^) 



(64) 



I 



d d q 1 



2tf 4 



d A q 



(27r)V(g 4 + 2tf 4 ) ' 



(65) 



The first term in expression (65) vanishes in dimensional regularization, while the second integral is 
convergent in the ultraviolet region by power counting in d = 4. The same argument is straightforwardly 
applied to the other Feynman diagrams which contribute at one-loop order to the the two-point 1PI 
Green's function of the gauge field A^. Thus, one can state that no divergent counterterms of the form 
"& 2 A'tJY^ show up at one-loop order, in agreement with the general expression of the counterterm given 
in eq.(54). Moreover, a similar analysis can be performed for the two-point mixed A-U Green's function 



in order to prove that no divergent counterterms of the type ^A^U^ are needed. In fact, the relevant 
one-loop Feynman integral in this case is of the type 



I(k) 



d d q -iV2$ 2 



(66) 



(2vr) d g 4 + 2i? 4 (q - k) 4 + 2tf 4 ' 
which is ultraviolet convergent by power counting. It is not difficult now, by making use of the decompo- 



sition (64), to extend these considerations to higher loop orders, so as to conclude that no counterterms 
proportional to the mass parameter "Q 2 are in fact needed in order to renormalize the theory, as expressed 



by equation (62 ). 



3.6 Renormalizability of the refined version of the replica model 



As already mentioned, the replica model, eq.( 11 ), admits a refined version pQ yielding a gluon propagator 
which does not vanish at the origin in momentum space, eq. The refined version of the replica model 
is obtained by adding to the action S rep n ca of eq.(ll) the following mass term 

,2 



rn 



d*x—(A«A; + U«Up 



which turns out to be BRST invariant on-shell, i.e. 



■"> I , I / n&S Te f 5S re f 

sS m = m z I d x ( c „ J + uj 



5b a 



Sb a 



where S re f stands for the refined action 



5, 



ref 



s, 



replica 



The mass term, eq.(67), modifies the behavior of the propagators according to 

k 2 + m? 



{A^k)A b u (-k)) 
(U«(k)U b u (-k)) 
(A a Jk)U b u (-k)) 



tab I e 

(k 2 + m 2 ) 2 + 2# 4 V 
k 2 + m 2 fe / 
(k 2 + m 2 ) 2 + 2tf 4 V 



iV2?? 2 



(k 2 + m 2 ) 2 + 2tf 4 



5 ab {5, 



i- uj 



kfj,k u 
k^ky 
k^k v 



(67) 

(68) 
(69) 



(70) 
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We notice that S m is left invariant by the discrete mirror symmetry ( 15 ) as well as by the rigid symmetries 
of eqs.(41) and (42). As the replica model, also the refined action S re f turns out to be renormalizable. 
To that purpose, we follow the previous set up and embed the mass term S m into an extended action by 
means of a BRST doublet of external sources (A, a) 



sA = a , sa = . 

It is thus easily verified that the expression S m is recovered from the BRST invariant term 

Sa = s J d 4 x [\{a%a% + u;u;) - | Xa) 

= I d 4 x (\a(A a ^ + Uffi) + \{Al d,c a + U; d,u a ) + | a 2 
when the sources (A, a) attain the values 



(71) 



Iphys" 



m 



A 



phys" 



0. 



(72) 



(73) 



The last term of the r.h.s of expression (72) is allowed by power counting, with £ being a constant 
parameter. For the complete extended action accounting for the mass term Sm, we write 

£' = £ + S a , (74) 



with S given in eq.(36). Besides the Ward identities (37)-(49), the extended action £' enjoys an additional 
symmetry 



d x — ic~— iuj"-=^- 

\ 5\ 5b a 5b a 

For the most general invariant counterterm we have now 



0. 



T,' = T, c + / d x 



(75) 



(76) 



where S c is given in expression (54) and where 04 stands for a free coefficient. This counterterm can be 
reabsorbed in the action £' by redefining the sources (a, A) and the parameter £ as 

Ao 



a = Z g Z A a, 



Z l g l 2 Z A l,i A , £ = [1 + e (2a + 2ai + 04)] £ , 



(77) 



in addition to the previous redefinitions (|58j) and ( 59 ) . This completes the proof of the renormalizability 
of the refined version of the replica model. 



4 Construction of BRST invariant local composite operators with the 
quantum numbers J PC = ++ ,2 ++ ,0 + and their z-particles content 

Glueballs are colorless composite gluon states which are classified according to the values of the angular 
momentum J, parity P and charge conjugation C. Within a quantum field theory framework, glueball 
states are constructed by means of suitable gauge invariant local composite operators carrying the quan- 
tum numbers J PC [9l [THl, [19] . The physical properties of the glueballs, i.e. masses, decay properties, etc., 
are thus encoded in the analytic properties of the correlation functions of the corresponding composite 
operators. 
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In what follows we shall look at the construction of BRST local invariant composite operators describing 
the lightest glueball states, J PC = ++ ,2 ++ ,0 h , and whose correlation functions do have good ana- 



lyticity properties. As already mentioned, see eqs.(26) and (27), this requirement can be achieved by 



looking at local operators which display an i-particles content. For the states with the quantum numbers 
J PC = ++ , ^ these operators are easily identified pQ and given by 

Oo++(x) = kwcw-^w^w). ( 78 ) 



0o-+ (a) = ftWW " U; u {x)U a pu {x)) . (79) 



The case of the state 2 ++ requires more care. Usually, the 2 ++ state is associated with the energy- 
momentum tensor which, in the present context, does not display the required i-particles content. Never- 
theless, we can construct a local BRST invariant operator which creates a pure 2 ++ state by demanding 
that it is symmetric, traceless and conserved. It has to be a symmetric tensor in order to have J = 2. It 
has to be traceless, otherwise its trace would give rise to a state with the quantum numbers of the scalar 
glueball ++ . Finally, it has to be conserved due to the fact that its divergence would be associated with 
a vector glueball state. To construct such an operator we follow the procedure already employed in |19j . 
and we consider the local invariant operator 

[02++WU = ( P ^ P vP ~ \ P ^ P «p) [K,i*)F$M ~ K a (x)U^(x)) , (80) 



where P pu = 5 pu d 2 — d p d u is the transverse projector. Expression (80) has the right properties for a pure 



2 ++ state. It is in fact symmetric, traceless and conserved. Let us also show that, as required, all three 



operators, eqs.(78), (79), (80), have an i-particles content. 
Introducing the i-particles field variables, eqs.([20|, 

K =7s( A m + <) ( 81 ) 

K =72^- O > ( 82 ) 

it is easily verified that 

O ++(x) = A" (ar)^" (a;) + higher order terms, (83) 

[02++O*OU = i P m P vp ~ \ P ^P a ^j (A^(x)^ CT (x) + < ff (x)A^(x)) + higher order terms, (84) 

O -+{x) = ^e^ pa (X a pu (x)r] a pa (x) + 'nf lu (x)X a pa (x)) + higher order terms, (85) 
where (A" , r?" ) stand for the quantities 

K» = d lx \?,-d v \% (86) 

= «-^<, (87) 

and where higher order terms in the fields have been neglected as they will not enter the evaluation of 
the correlation functions at one-loop. 
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4.1 Obtaining the Kallen-Lehmann spectral representation for the one-loop corre- 
lation functions 



At one-loop order, the two-point correlation functions of the glueball operators, eqs.(83), (84), (85), take 
the general form 



(Oi(k)Oi(-k)) 



d 4 p 



1 



1 



;fi(p,k-p) ■ 



(88) 



i-ioop J (2vr) 4 (k-p) 2 + iV2$ 2 P 2 -iV2$ 2 ' 
where fi(p,k — p), i = ++ ,2 ++ ,0 , are polynomials in the scalar products of the momenta (k,p). 



There are several ways to show that expressions of the form ( 88 ) do exhibit in fact a Kallen-Lehmann 
spectral representation, namely 



(Oi(k)Oi(-k)) 



l—loop 



T + k 2 



(89) 



One possibility is to follow the procedure employed in [6] , relying on the use of Feynman parameters and 
of dimensional regularization. Here, we shall follow another route and establish the spectral representa- 
tion by making use of the Cutkosky's rules in Minkowski space. At the end we shall perform an analytic 
continuation to complex masses in Euclidean space. We shall explicitly see that the results obtained this 
way turn out to coincide in fact with those of [6]. 

Cutkosky's rules enable us to obtain the discontinuity of Feynman amplituted across the cuts. For 
that, each Feynman propagator corresponding to a cut internal line is put on-shell, according to 



1 



p 2 — m 2 



2tt9{p )8(p 2 - m 2 



Thus, for the discontinuity of the one-loop Feynman integral in Minkowski space 

d 4 p 1 1 



(2tt) 4 (k-p) 2 -m 2 p 2 



rrio 



Ji(p,k-p) , 



we get 



ImTi 



d A p 



(27r) 2 9((k - p)°)6((k - p) 2 - m 2 )6(p°)5(p 2 - m 2 )f t (p, (k - p)) . 



(90) 



(91) 



(92) 



We are ultimately interested in the case where the masses involved are complex, i.e. m 2 — > m\ = i\[2 r Q 2 
—i^/2'd 2 with the momenta defined in Euclidean space. We then evaluate expression 
(92) as a function of real masses and Minkowsky momenta and we analytically continue the resulting 



and m 2 



rat 



expression to complex masses in Euclidean momentum space. This leads to the following prescription for 



the expression of the spectral functions entering expression ( 89 ) 



p r : ep (T,m x ,m v )=pt Imk (T,m u m 2 ) 



where pf {t, m\, m?) in the right hand side is obtained from expression (92), namely 

1 



P 



Mink 



TT 



(93) 



(94) 
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The integral (92) is readily evaluated, and we find the spectral representation at one-loop order for 
each of the two-point correlation function of the glueball operators: 



(O 0++ (k)O 0++ (-k)) 



1—loop 



2{N 2 - 1} 
8^2 



dr 



1 



¥/8^ 2 

T + k 2 1 2 +T 



(95) 



<[02++(*)U [<V+(-*)U 



1—loop 



\(N 2 



8ir 2 



1) r°° \/l 

hV2^ r + k 2 



^(8^ 4 ) 2 r 2 + 2(8^)T 4 + ^ 6 
8 2 



(96) 



(O - + (k)O - + (-k)) 



1—loop 



8(N 2 - 1) 
8^2 



x \l 1 

1 t + k 2 



t2 



8tf 4 ) 



(97) 



where the threshold is found to be given by the expression 

Notice that all spectral densities p[ ep , i = ++ ,2 ++ ,0 h are positive in the range of integration. It is 



also worth pointing out that expressions (95) and (97) coincide with those already reported in [6]. 



5 Establishing a SVZ-type sum rules. An exercise on the spectrum of 
the glueballs 

Having identified a good set of BRST invariant local composite operators with quantum numbers J PC = 
++ ,2 ++ ,0 ^ and with positive spectral densities, we might attempt at achieving a qualitative prelim- 
inary analysis of the ratios of the three masses m 2 ++ , m 2 , ++ , "2-Q-+ , in order to obtain a first indication 
of their location. To that purpose, we shall work out a kind of phenomenological SVZ-type sum rules 
containing a free coefficient a which parametrizes our lack of knowledge of the true physical spectral 
functions. Further, an expression for the masses of the three glueballs ++ ,2 ++ ,0 ^ as a function of 
the coefficient a is obtained by making use of the Borel transformation. In the spirit of the SVZ sum 
rules, we shall look thus at the existence of an interval for the parameter a, for which the location of the 
three masses m 2 ]++ ,m 2 t ++ ,m 2 ] _ + is in agreement with the lattice data, see ref.[9], according to which the 
lightest state is the ++ , followed by the 2 ++ , the ^ being the heaviest one. The existence of such an 
interval will be taken as a first encouraging evidence for a more complete quantitative analysis. Let us 
thus proceed by describing how our phenomenological SVZ-type sum rules are worked out. 

As it is customary in the SVZ approach to QCD [10], we start by considering the two-point correla- 
tion functions 

Mi) = J d 4 xe^(O i (x)O i (0)) , (98) 

where Oi, i = ++ ,2 ++ ,0 , stand for the local composite gauge invariant operators which generate 
glueball states with quantum numbers J PC = ++ , 2 ++ , ^ . 

On physical grounds, a truly nonperturbative evaluation of Hi(q 2 ) would enable us to write an exact 
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Kallen-Lehmann spectral representation 



2 . , r . ImUdr) 



U i (q I ) = - dr V , (99) 
TT ,/n r + q z 



which is expected to follow from the unitarity and analyticity properties of the underlying nonperturbative 
theor>|^] We thus proceed by employing a one-resonance parametrization for Imllj(r) [2D], namely 

= n . 6{T _ m |) + fl(r _ T * )p f y- (T) (100 ) 

where m| denotes the glueball mass in the i-th channel, Tq the threshold for the continuum part of the 
spectrum and f^ hys (T) the corresponding positive spectral density. Of course, the real values of IZi, mf, 
Tq and of the spectral density p^ hys {T) are unknown. So far, the best which can be done is computing the 



correlation functions (98) by trying to encode as much nonperturbative effects as possible, as summarized 



by the following equation 



Ui + /°°dr^M = n7, (ioi) 



q 2 + m 2 L T + q 2 



where II stands for the expression of the correlation function ( 98 ) which one has been able to evaluate 



in practice. Expression (101) establishes the so-called sum rules, enabling us to give estimates of the 
glueball masses in terms of the nonperturbative parameters present in IL P . 

As shown in their seminal work [T2l [T3l 114] , this approach has been proven to be rather success- 
ful in order to obtain estimates of hadronic masses in QCD. 

Here, we shall not follow the standard path of the SVZ sum rules. We remind here that in the SVZ 
sum rules approach the quantity II™ P in the right hand side of eq.(lOl) is obtained through several steps. 



First, one starts from the usual QCD Lagrangiaij 10 | which is employed to evaluate the perturbative con- 



tributions till a certain order. Further, non-perturbative contributions are included by the introduction 
of a suitable set of condensates obtained via operator product expansion (OPE). Finally, additional non- 
perturbative effects are obtained from instantons, see, for example, |20l 121] for a detailed account. In the 
present work, we follow a slightly different route. The non-perturbative effects are already encoded in 
the starting action by means of a non-perturbative mass parameter t? 2 which accounts for the presence 
of the Gribov horizon. More precisely, we shall attempt attempt at evaluating the right hand side of 



eq.(lOl) by employing a trial nonperturbative action S np fulfilling the following requirements: 



i) the action S np exhibits the nonperturbative effects of the Gribov horizon, whose presence is 
parametrized by the Gribov mass parameter, here denoted by i? 2 . Moreover, in the limit $ 2 = 0, 
in which the Gribov horizon is removed, the action S np ($ 2 ) reduces to the usual perturbative 
Faddeev-Popov action Sfp, i-e. 

S np ({) 2 ) = S FP , (102) 
$ 2 =o 

ii) it accounts for gluon confinement. This means that the two-point correlation function of the 
elementary gluon field appearing in S np ($ 2 ) cannot be cast in the form of a spectral representation 
with positive spectral density, so that it cannot be interpreted as the propagator of a physical 
particle. 



We remind here that, in some cases, the spectral representation, eq.(98l, might require appropriate subtraction terms 
in order to ensure convergence. These terms are not written down, as they will be removed once the Borel transformation 
will be taken. 

w i.e. J d A x {\F 2 + V^mA^ + mi>lp). 
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• Hi) it is renormalizable, meaning that consistent calculations can be worked out. 



iv) it enables us to introduce gauge invariant or, equivalently, BRST invariant local composite oper- 
ators 0™ p with the quantum numbers J FC = ++ , 2 ++ , whose two-point correlation functions 



(103) 



exhibit the Kallen-Lehmann spectral representation, a feature which we shall consider here at one- 
loop order only, i.e. 

n"V,tf 2 ) = / _ dr PtSLL + (h 2 ) . (104) 



Requirements 



-iv 



look quite stringent. As such, they might provide a satisfactory set up in order to 



achieve a useful predictive expression for the correlation functions II™ p (q 2 ,'!? 2 ) 



eq.(103). 



Two examples of field theories compatible with the requirements i)—iv) are known so far. The first 
example is provided by the Gribov-Zwanziger action [21 O H] and its refined version |5J. Though, till 
now, the item iv) has not yet been completely settled [6]. The second model available is the replica 
model discussed here. As already mentioned, this model can be considered as a useful renormalizable 
laboratory reproducing most of the features of the Gribov-Zwanziger theory. It contains a soft mass 
parameter r & 2 which plays a role akin to that of the Gribov mass parameter, while allowing for the in- 
troduction of BRST invariant local composite operators with quantum numbers J = ++ ,2 ++ ,0 
whose two point correlation functions exhibit the Kallen-Lehmann spectral representation at one-loop 
order, thus fulfilling item iv). It is worth mentioning that both models rely on the mechanism of the soft 
breaking of the BRST symmetry. In the following, the action of the replica model will be taken as our 
trial nonperturbative action, i.e. 

S ''(l? ) — S re pH ca . (105) 



According to expression (101), we shall thus elaborate on the equation 



q 2 + mf 



+ 



oo n phys 

dr 



T + q 2 



2 V / 2# 2 



n re P( T ) 

dr P -^l + <D(h 2 ) 
T + q 2 



(106) 



where p^ v {r) are the one- loop spectral densities given in eqs.(95),(96),(97). 



As in the SVZ sum rules approach [10| . before going any further we need to provide an estimate for 
the quantity 



dT 



pf VS {r) 
T + q 2 



appearing in eq.(106). Concerning the threshold Tq, we shall write 

4, = a2\f2-d 2 , i = 0++, 2 ++ , 0" 



(107) 



(108) 



where a > 1 stands for a free parameter accounting for the difference between the value of the physical 
threshold, Tg, and that which we have been able to evaluate in practice, i.e.: 2^f 7 2 r & 2 . One should notice 
that, in principle, the physical threshold Tq could be different for each glueball correlation function. 
Instead, we are employing a unique parameter a for all states i = ++ ,2 ++ ,0 ^, a choice which is 



motivated by the fact that three correlation functions (95), (96), (97) exhibit the same threshold. Further, 
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for what concerns the spectral density p? hys (r), we shall proceed by assuming that, at the one-loop order 
approximation considered here, it can be well approximated by p^ (t), and that possible deviations start 
to show up at higher orders, i.e. we shall set 



pt VS (r) = ft> (r) + 0{tf) 



Therefore, for the final form of the sum rules, we write 

r2V2atf 2 



q 2 + m\ 



rep I 



d.T 



2s/2$ 2 



r + q 2 



0(k 2 



(109) 



(110) 



Equation (110) will be employed for a qualitative discussion of the ratios of the masses of the lightest 



glueball. This will be done by looking at the existence of an interval for the phenomenological parameter 



a for which the ratios of the glueball masses obtained through eq.(llO) are in qualitative agreement with 



the available lattice data. Let us also underline that the necessary non-perturbative ingredient for a 
study of the spectrum of the glueballs, even at the qualitative level, is encoded here in the presence of 
the parameter i9 2 which, to some extent, accounts for the effects of the Gribov horizon. 



5.1 A formula for the glueball masses 



Expression (110) enables us to extract a formula for the glueball masses. To that purpose, we shall first 

(111) 



compute the Borel transformation of the first moment of the sum rules (110), namely 



B 



M 



q^ + m, 

where Bu stands for the Borel operator [9"1 [TO] 
BAi(f(q 2 )) = hm ni9 




2V2 a-d 2 



dr 



p7 p (t) 



T + 



4)n 



111 



(9 



2\n+l 



1 fd n f(q 2 



d n q 2 



=M 2 =fixed 



and M 2 is the so called Borel mass [H [10] . 
From 



B 



M 



q 2 + mj 



e m 



one easily gets 



c\ i 

TZjiri; e m 2 



2y n 2a-d 2 



2s/2& 2 



dr t p i p (r) e 



Taking now the Borel transformation of eq.(110) 



2v^2 ad 2 
2V2i? 2 



drp[ ep (r)e ~m 1 



(112) 



(113) 
(114) 



(115) 



and evaluating the ratio between eq.( 114[ ) and eq.(115), we get the mass formula we are looking for 

r2V2ai9 2 



VTLa 



2 v / 2)? 2 



drr p- (r) e 



r 2 v / 2aiS 2 , rep, x 
hV20 2 dT Pi ( T ) e ^ 



(116) 



An interesting aspect of expression (116) is that the mass mf exhibits an explicit dependence from both 



threshold 2^/2-d 2 and spectral density p r ^ v (j} 
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5.2 A qualitative analysis of the ratios of the glueball masses rrii 



Q++i m 2++i " b o-+ 



We have now all ingredients for a qualitative analysis of the glueball masses m^ ++ , m 2 4 



use of expressions (|95|),(|96|),(|97|), equation (j 1 1 6|) can be written in explicit form giving 

8tf 4 



m 



o++ 



dr-y/r 2 " 



8t? 4 



+ r z e m2 



2^2 ai? 2 
2V2i? 2 



dr- 



e 



. . Making 



(117) 



m 



2++ 



2 v / 2ai? 2 



8t? 4 (f (8tf 4 ) 2 r 2 + 2(8t? 4 )r 4 + §r 6 ) e" 



r2\/2ai9 2 , 



-8i? 4 (§(8i? 4 ) 2 t 2 +2(8i? 4 )t 4 + § r 6 ) 



(118) 



m 2 _+ 



J2 2 0f/^P-^ 4 ) 3/2 

f 2v^a^ 2 j (r 2 -8i? 4 ) 3/2 
J2v / 2^ 2 aT 



e m 2 



(119) 



It is possible to rewrite these expressions in such a way to make all the masses proportional to the 
mass parameter i? 2 . We introduce the variable t = 2v ^ 2 an d the parameter p = 2 ^Qt ■ The resulting 
expressions read 



m 



0+H 



(a,p) 



-pt 



(120) 



m 



2++ 



[a,p) 



tfdtV? 



ia t 2 + 2t 4 + 3 t 6 )e - pt 



-Pt 



(121) 



Iidt' 



3 / 2 e - P t 



^, (* 2 -l) 3/2 c - P t 



(122) 



The masses become now functions of the parameters a and p. The output of our results are shown 
in Fig.l, Fig. 2 and Fig. 3. From Fig.l one can see that, when the parameter a belongs to the interval 
1 < a < 1.8, the masses of the three lightest states are in qualitative agreement with the available data, 
i.e. m,Q ++ < m 2 ++ < mg_ + , a feature which holds for all values of p, as shown in Fig. 2. and in the three 
dimensional plots of Fig. 3. The existence of the interval 1 < a < 1.8 is regarded as a preliminary test 
towards a more quantitative analysis of the glueballs spectrum. 
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Figure 2: Glueball masses as functions of p for a = 1.3 
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Figure 3: Glueball mass ratios 



6 Conclusion 



In this work several aspects of the replica model introduced in [TJ have been addressed. The relevance of 
this model relies on the fact that it reproduces several features of the Gribov-Zwanziger theory [21 El S] . 
As such, it can be employed as a useful laboratory to investigate, within the framework of the Euclidean 
quantum field theory, the consequences of having a confining gluon propagator of the Gribov type, eq.Q. 

We have shown that the replica model is a renormalizable theory, while enabling us to introduce suitable 
BRST invariant local composite operators with the quantum numbers of the lightest glueball states, 
J PC = ++ ,2 ++ ,0 '". The corresponding correlation functions have been evaluated at one- loop order 
and shown to display a spectral representation with positive spectral densities. 

A first check of the fact that the nonperturbative effects related to the Gribov horizon might be useful for 
the investigation of the properties of the glueball spectrum has been achieved through a phenomenological 



SVZ-type sum rules, eq.(HO). Although the analysis presented here is only at the qualitative level, it can 
open the road for a more quantitative study in which the most recent available lattice data on the gluon 
two-point correlation function can be employed to get numerical estimates of the Gribov mass parameter 
and of the related dimension two condensates |22J. This could enable us to obtain quantitative estimates 
for the glueball masses Tn^ ++ , m 2 , ++ , m 2 ] _ + , to be compared with the available data. 

Finally, although not reported, we mention that plots similar to those of Figures 1,2,3 can be obtained 
in the case of the decoupling propagator, eq.Q, which also displays the z-particles decomposition 

k 2 + m 2 If 1 1 \ , . 

+ ,, . , ■ • (123) 



(k 2 + m 2 ) 2 + 2i? 4 2\k 2 + m 2 + i^ 2 k 2 + m 2 - iV2^ 
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